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We present calculations of spin-dependent electron transport through single organic molecules 
bridging pairs of iron nanocontacts. We predict the magnetoresistance of these systems to switch 
from positive to negative with increasing applied bias for both conducting and insulating molecules. 
This inverse magnetoresistance phenomenon does not depend on the presence of impurities and 
is unique to nanoscale magnetic junctions. Its physical origin is identified and its relevance to 
experiment and to potential technological applications is discussed. 
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I. INTRODUCTION 

Electronic devices composed of ferromagnetic elec- 
trodes separated by a connecting junction whose resis- 
tance can be altered with magnetic fields offer enormous 
potential for technological applications such as mag- 
netic random access- memory (MRAM). Spin-dependent 
transport in ferromagnetic nanosystems is currently re- 
ceiving increased attention owing to the expectation 
that nanoscale size constraints may maximize the po- 
tentially useful magnetoresistance effects. In particu- 
lar, theories have predicted that molecules bridging fer- 
romagnetic electrodes may present pathways for build- 
ing ultra-small spintronic devices that take advantage 
of weak molecular spin-orbit interaction to maximize 
the magnetoresistancei2i2i4 Recently, experiments have 
been reported in which electrons transmitted through in- 
sulating octanethiol molecular self-assembled monolayers 
(SAMs) between nickel contacts^, and through lOOnm 
organic molecular films with magnetic electrodes^'^, re- 
tained their spin polarization leading to observed magne- 
toresistance, proving that individual or groups of organic 
molecules may potentially be used as building blocks for 
the ultimate miniaturization of spintronic devices. How- 
ever, bias dependence and inversion of the magnetoresis- 
tance were also reported and were attributed^ to the pres- 
ence of impurities within the molecular junction. Bias 
dependence of the magnetoresistance has since been pre- 
dicted in insulating molecular junctions^ , indicating that 
bias-dependent magnetoresistance may be intrinsic to 
molecular spintronic devices. However, most technolog- 
ical device applications based on magnetoresistance re- 
quire a relatively bias-independent magnetoresistance for 
stable device operation. In this article, we report predic- 
tions of magnetoresistance in ideal molecular spintronic 
junctions where single molecules bridge iron electrodes, 
systems that until now have received little experimental 
or theoretical attention. We predict strong bias depen- 
dence in the positive magnetoresistance regime, but we 
also predict the magnetoresistance to become negative 
with increasing applied bias and we show that the mag- 
netoresistance is potentially larger and stable under ap- 
plied bias in this negative magnetoresistance regime. Our 
calculations reveal negative magnetoresistance to be an 



intrinsic property of some molecular junctions (i.e., un- 
like the effect reported in Ni/octanethiol/Ni junction£^«, 
not dependent on the presence of impurities) and thus 
offer a path to stable operation of devices (such as molec- 
ular MRAM) based on molecular spintronics but operat- 
ing in the negative magnetoresistance regime. Therefore, 
understanding the origins of negative magnetoresistance 
is paramount to fundamental studies of spin-dependent 
transport in molecular junctions and potentially applica- 
ble to technological endeavors. 

We explore theoretically two very different molecular 
junctions: a conducting (conjugated) benzene-dithiolate 
(BDT) molecule and an insulating octane-dithiolate 
(ODT) molecule bridging ferromagnetic Fe nanocontacts. 
The phenomena that we predict appear in both. 



II. THEORY 

In our calculations the system is partitioned into semi- 
infinite ideal ferromagnetic source and drain leads and an 
extended molecular junction consisting of the molecule 
and clusters of nearby Fe atoms^; for BDT this "ex- 
tended molecule" is shown in the inset of Fig^Ja). The 
Fe-thiol bonding geometries are estimated with ab initio 
relaxationsi2i& The electronic structure of the Fe clus- 
ters is described by a tight-binding Hamiltonian and 
non-orthogonal s,p,d basis. The tight-binding param- 
eters are based on fits to ab initio band structures of 
Fe crystals^-'' and have previously been employed suc- 
cessfully to study magnetic multilayers^^, magnetioi^ii^ 
and non-magnetioi^ nanocontacts, and ferromagnetic 
Ni molecular junctionsi. For Fe structures with sur- 
faces this model yields spin-resolved surface densities of 
states and enhanced surface magnetic momentaii sim- 
ilar to those of ab initio surface calculationsi^. Thus 
our model incorporates both the bulk and surface mag- 
netic properties of Fe which together influence the spin- 
dependent transport through a nanoscale junction of bulk 
Fe leads bridged with molecules. Spin flip processes are 
not considered, consistent with weak spin-orbit interac- 
tion in molecules and the high degree of spin polariza- 
tion retained in spin-dependent observations of molec- 
ular junctions .-5*^ The molecular electronic parameters 



are described by a tight-binding formalism based on 
extended-Hiickel (EH) theory. This approach has been 
used successfully to explain the experimental current- 
voltage characteristics of molecular nanowires connect- 
ing non-magnetic metal electrodes»iS*i^4S£ The EH pa- 
rameters are based on atomic ionization energies while 
the electronic parameters from Ref. describing the 
Fe clusters are defined up to an arbitrary additive con- 
stant. We adjust this constant to align the Fermi energy 
of the contacts relative to the highest occupied molec- 
ular orbital (HOMO) of the BDT molecule, according 
to the difference in ionization energies of the isolated 
structures estimated with ab initio total energy calcu- 
lations. For Fe electrodes and the molecules consid- 
ered here, this agrees well with the difference between 
the work function of Fe and the HOMO energy of the 
isolated molecule obtained from density functional the- 
ory, a method that has been used successfully to align 
the Fermi energy of gold with the HOMO of BDT.^^ 
The potential profile of the nanocontact arising from 
the bias voltage applied between the two electrodes is 
difficult to calculate from first principles since it is a 
non-equilibrium many-body property. However, appro- 
priate heuristic models for the profile can yield accu- 
rate results for the current We adopt this approach 
here, assuming the majority of the applied bias to drop 
over the metal-molecule interface, ^^'^■^•^^^^ Nearly identi- 
cal current- voltage characteristics emerge if the details of 
the model potential profile are varied. Thus these details 
are not crucial to our predictions. Our transport calcu- 
lations are based on Landauer theorjiS^ and Lippmann- 
Schwinger and Green's function techniques. We define 
the junction magnetoresistance as J MR — / . 

Ipar (lanti) IS the Current through the junction for paral- 
lel (anti-parallel) magnetizations on the Fe electrodes. 

III. RESULTS 

Fig. n^a) shows the calculated transmission probabil- 
ities at zero bias for the BDT molecular junction. The 
spin up — > up (uu) and down — > down (dd) transmis- 
sion is for parallel magnetization of the two contacts; 
up —^ down (ud) is for anti-parallel magnetization. Quite 
similar results for the transmissions and currents are ob- 
tained if the 5x5 outer Fe layer is deleted from each of 
the clusters indicating that the clusters and ideal leads 
in our model adequately represent real (macroscopic) Fe 
leads. The molecular HOMO and HOMO-f give rise 
to the strong transmission 1.5-2eV below the Fermi en- 
ergy Ep = OeV; the peaks 3-4eV above the Fermi en- 
ergy are due to the molecular LUMO and LUMO-|-f. 
The sharp resonance at 3.1eV results from extremely 
weak hybridization between the Fe leads and the iso- 
lated molecule's LUMO which has very little sulfur con- 
tent. The transmission resonance is present and trans- 
mits for all spin configurations. More important, how- 
ever, for potential device applications are resonant states 
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FIG. 1: (Color online) (a) Transmission probabilities vs. en- 
ergy at zero bias for the spin up — > up (uu), down down 
(dd), and up —> down (ud) configurations of the BDT molec- 
ular junction shown in the inset. Fermi energy = OeV. (b) 
Current vs. voltage for parallel (p) and anti-parallel (a) mag- 
netizations. JMR (J) is also shown. 



arising from strong hybridization between the molecu- 
lar sulfur and the Fe d-electron orbitals. In Fig. QJa) 
these states give rise to the broad transmission features 
between the Fermi energy and the HOMO resonances 
mentioned above. They allow electrons to transmit effi- 
ciently through the metal-molecule interface and so give 
rise to moderately strong transmission dominating the 
current in the experimentally accessible moderate bias 
regime. Together with the energetically localized and 
highly spin-split Fe d bands, these hybrid states give rise 
to strongly spin-dependent transport characteristics. 

At zero applied bias for parallel magnetization of the 
contacts, the spin up Fe d orbitals hybridize very strongly 
with sulfur states on both sides of the molecule as shown 
for a typical state in Fig. |2Ia). This results in the 
dominant uu transmission near (below) the Fermi en- 
ergy in Fig. n^a). The spin down d orbitals, occurring 
higher in energy, hybridize only weakly with molecu- 
lar sulfur states. For anti-parallel magnetizations [i.e., 
up{down) — > down[up) transmission], only the source 
(drain) contact hybridizes strongly with the molecular 
end group (Fig. [2Ib)), resulting in much weaker ud trans- 
mission (anti-parallel magnetization) than the uu trans- 
mission in Fig. n^a) . As the hybrid states are not located 
exactly at the Fermi energy they contribute to the current 
at very low bias through their off-resonant tails. 

An applied bias breaks the symmetry between the 
source and drain electrodes in the parallel magnetization 
configuration. Thus the symmetric zero bias molecular 
orbitals (as in Fig. [2Ia)) are supplanted by states associ- 
ated mainly with either the source or drain electrode as 
shown in in Figure [2Ic,d) for the uu transmission. The 
orbital in Figure[2Ic) is shifted upwards in energy and is 
associated with the source electrode, while that in Fig- 



3 



(a) 



(c) 



(uu OV) (ud OV) 

(uu Source 0.1V) (uu Drain 0.1 V) 



(e) 




(f) 



(uu 1.0V) 




(ud 1 .OV) 



FIG. 2: (Color online) Representative hybridized energy 
eigenstates of the entire extended molecule near the Fermi 
energy are depicted on the BDT molecule and the 2x2 Fe 
planes on either side, (a) At OV the uu source and drain elec- 
trode d states hybridize with sulfur orbitals on both sides of 
the molecule forming transmitting states, (b) Only the source 
strongly hybridizes with the sulfur states in the ud configu- 
ration. At O.IV, uu states in (a) split into a state mainly 
on the source (c) that shifts upwards in energy and a state 
mainly on the drain (d) that shifts downwards, (e) At l.OV 
uu hybridized states on the source and drain have shifted well 
apart, (f) Source up hybridized states are shifted into reso- 
nance with drain down states in the ud configuration creating 
states that transmit strongly through the molecule. 



ure|2Id) is shifted downwards and is associated with the 
drain. These orbitals can be traced to the zero bias or- 
bital in Figure[2Ia) and their energy splitting corresponds 
to the applied bias. Since the energies of these states are 
not perfectly aligned, simultaneous resonant transmission 
through both ends of the molecule is impossible. How- 
ever, since the states remain close in energy at moderate 
bias, they both contribute coherently to near-resonant 
transmission via the extended molecule's Green's func- 
tion in the Lippmann-Schwinger scattering formalism: 
There is appreciable coherent mixing between their off- 
resonant transmission tails resulting in relatively strong 
uu transmission. The corresponding molecular orbital for 
ud transmission resembles that shown for the source elec- 
trode in the uu configuration in Figure \^c) , however a 
state corresponding to the uu drain state is absent. The 
reverse is true for du transmission and since these two 
transmission processes are mutually incoherent, there is 
no mixing of their off-resonant tails and the current in 
the case of anti-parallel magnetizations is weaker. 

Fig. ni^b) shows the corresponding finite-bias current- 
voltage characteristics. At low bias, the current grows 



as more hybridized states contribute to the transmission. 
Since the uu transmission dominates, the parallel magne- 
tization current is larger, resulting in a positive magne- 
toresistance of about 10 percent. At larger applied bias, 
the applied electrostatic potential breaks the symmetry 
between the contacts further and the uu transmission de- 
creases as the hybridized Fe d states associated with the 
source and drain move apart. However, at high bias in 
the case of anti-parallel magnetizations, hybridized states 
associated with down d electron orbitals on the drain 
are shifted into resonance with up d electron states on 
the source and hybridized transmitting states connecting 
both sides of the molecule with the electrodes form in the 
case of up down transmission as shown in Fig. Olf), 
enhancing the transmission. No such resonant effect oc- 
curs for du transmission which decreases with increasing 
bias and so the overall effect of the applied bias on the net 
transmission for anti-parallel magnetization is small. No 
resonant enhancement occurs at increased bias for par- 
allel magnetizations where hybridized states are shifted 
apart (Fig. [2Ie)) and transmission is reduced. Therefore, 
because it selectively depresses transmission for parallel 
magnetization of the contacts, the bias-induced symme- 
try breaking results in a crossover with increasing bias 
from positive to negative values of the JMR in Fig. db). 

The wave functions of the hybridized states responsi- 
ble for transport at moderate bias are located primarily 
on the Fe electrodes. Therefore their energy levels are 
pinned to the electrodes and shift rigidly with applied 
bias explaining why our moderate bias results are insen- 
sitive to the assumed potential profile. If one partitions 
the entire system into halves associated with either the 
source or with the drain so that the hybridized states on 
the source (drain) side of the system are grouped with 
the source (drain) to form a "composite" source (drain) 
electrode, the negative magnetoresistance argument can 
be rephrased in the language of conventional solid-state 
magnetic tunnel junctions: At the bias required to bring 
hybridized down states on the drain into resonance with 
hybridized up states on the source, the spin polarization 
of the composite drain is negative (the current is dom- 
inated by down electrons) while that of the composite 
source remains positive. Therefore, the JuUiere^^ theory 
of tunneling would predict negative magnetoresistance. 

Because its mechanism relies on the simultaneous sup- 
pression of transmission for parallel magnetizations on 
the contacts and enhancement of the transmission in 
the anti-parallel ud transmission channel, this intrinsic 
negative magnetoresistance requires strongly spin-split 
energy bands providing a large magnetic moment. No 
such negative magnetoresistance has been reported for 
calculations on molecules bridging ferromagnetic nickel 
junctions^'^''^'^ and our calculations on similar systems 
suggest that this different behavior may be related to the 
smaller spin-splitting in Ni: At sufficient bias to induce 
resonant transmission in the ud channel (the required 
bias is relatively small for Ni) much symmetry and over- 
lap remain in the uu transmission channel where energy 



levels are not adequately shifted apart. Therefore the 
mechanism should be less effective in nanosystems based 
on nickelSSSa. 

For BDT bridging Fe nanocontacts, the crossover to 
negative magnetoresistance occurs at relatively low bias, 
near 0.7V, well within the range of current experimental 
techniques and therefore the predicted negative magne- 
toresistance should be accessible to experiment. With 
magnitude reaching a value of more than 30 percent, the 
negative magnetoresistance may be more easily observed 
in this system than positive values. Clearly the mecha- 
nism is expected to be unique to systems that can both 
sustain a significant potential drop and allow for reso- 
nant transmission of d-electrons through the junction, 
i.e., nanosystems. Negative magnetoresistance emerging 
as a result of applied bias has already been observed ex- 
perimentally in systems where ferromagnetic electrodes 
are separated by thin film tunnel junctions of only a few 
nanometers in width where the junction couples strongly 
to the electrode d states^LS. Since d-electrons are not 
expected to contribute to transport in conventional thin 
film magnetic tunnel junctionsiSi^iMiS^, negative magne- 
toresistance when observed in those devices is attributed 
to impuritieS]i2& However, we predict intrinsic negative 
magnetoresistance in ideal systems demonstrating that 
bias-dependent and negative magnetoresistance are not 
necessarily signatures of impurities in molecular nanosys- 
tems. 

Interestingly, the region exhibiting the greatest 
bias independence in our magnetoresistance calcula- 
tions emerges in the negative magnetoresistance regime 
(around 1.7V in Figure QJb)). Since magnetoresistance 
devices can in principle exploit magnetoresistance of ei- 
ther sign but require a relatively bias independent mag- 
netoresistance for stable device operation, molecular-Fe 
systems may be attractive candidates for technological 
applications in devices based on molecular spintronics 
(molecular MRAM) but operated in the negative mag- 
netoresistance regime. 

For comparison, we have also studied an Fe/BDT/Fe 
molecular junction having a different geometry where sin- 
gle Fe tip atoms have been added to each of the clusters 
and the molecular sulfur atoms bind directly over them 
at the same heights considered previously. This geome- 
try results in stronger and sharper resonances associated 
with the hybridized states and therefore the mechanism 
responsible for negative JMR is even more prominent. 
Here we find the cross-over to negative magnetoresistance 
to occur at a lower bias voltage of Q.25V. 

Next we consider ferromagnetic Fe contacts bridged 
with an octane-dithiolate (ODT) molecule, shown in Fig. 
|3fa) with the zero bias transmission probabilities. As 
ODT is an insulator, the strongest transmission features 
near the Fermi energy are weak resonant features asso- 
ciated with the HOMO, about 1.8eV below the Fermi 
energy. Other weakly transmitting states appear closer 
to the Fermi energy (about 0.5 eV below) due to hy- 
bridization of Fe d and molecular sulfur orbitals. These 
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FIG. 3: (Color online) Calculated spin-dependent transmis- 
sion, currents and JMR of the ODT molecular junction shown 
in the inset. Notation as in FiglTI 



states are similar to those predicted in the case of BDT: 
the up — > up transmission is largest as hybridization oc- 
curs at both ends of the molecule allowing for efficient 
transmission into and out of the metal-molecule interface, 
decreasing the width of the tunnel gap. The calculated 
current- voltage characteristics are shown in Figure ISfb). 
At moderate bias, the symmetry between the source and 
drain electrodes in the parallel magnetization configura- 
tion leads to positive magnetoresistance. At higher ap- 
plied bias, by the same mechanism described for BDT, 
the magnetoresistance becomes negative, but here the 
crossover doesn't occur until nearly 1.2V. Again a large 
region of relative bias stability emerges in the negative 
bias regime (around 1.5V in Figure Ol^b)) with possible 
implications for future molecular spintronic technology. 



IV. CONCLUSION 

In conclusion: We have predicted an inverse magne- 
toresistance phenomenon intrinsic to systems where fer- 
romagnetic Fe electrodes are bridged by two very dif- 
ferent molecules (conducting and insulating), and for 
different bonding geometries. We also predict qualita- 
tively similar behavior in Fe junctions bridged by a pair 
of Fe atomsii^ Thus negative magnetoresistance should 
be available in a variety of molecular junctions under 
appropriate conditions. The value of the bias required 
to induce negative magnetoresistance depends on the 
chemistry of the junction: G.25-0.7V for Fe/benzene- 
dithiolate/Fe (depending on the bonding geometry), 1.2V 
for Fe/octane-dithiolate/Fe and higher stilU^ for Fe con- 
tacts bridged by a Fe dimer. These changes in the chem- 
istry make it possible to vary the onset of negative mag- 
netoresistance, the magnitudes of both positive and neg- 
ative magnetoresistances and the regions of bias indepen- 



5 



dence important for device applications. 
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